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We analyse pionic couplings of heavy mesons combining PCAC with the dispersion quark model 
to calculate the relevant transition form factors. Ground states and radial excitations are considered. 
For the ground state coupling constants the values g — 0.5 ± 0.02 in the heavy quark limit, and 
gB*B-n = 40 ± 3, go'D-n = 16 ± 2 are obtained. A sizeable suppression of the coupling constants 
describing the pionic decays of the radial excitations is observed. 
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Pionic coupling constants of heavy mesons are basic quantities to understand the heavy meson lifetimes: they 
govern strong decays of heavy mesons by the emission of a pion which in practice turn out to be most important 
0^ ■ strong decay modes. The knowledge of the pionic coupling constants of the ground-state as well as the excited heavy 
mesons is timely in particular in view of the new experimental data on the excited heavy mesons |l]] . 

The heavy quark (HQ) symmetry provides important relations between the coupling constants of pionic transitions 
of various mesons containing a heavy quark. In the leading I/tuq order (LO) the heavy quark spin decouples from 
other degrees of freedom Q. Thus in the limit of an infinitely heavy spectator quark, the respective pionic decay 
rates for mesons with different spins but same quantum numbers of the light degrees of freedom (total momentum of 
. the light degrees of freedom j, angular momentum L, etc) || are equal. 

Note, however that the heavy quark expansion which provides many rigorous results for semileptonic transitions 
between heavy mesons, turns out to be much less efficient in the case of the strong transitions. In fact, these two types 
H ' °f processes have quite different dynamics. Semileptonic transitions of heavy mesons are induced by the heavy quark 
^ | weak transition and in this case HQET predicts not only the structure of the 1 jrriQ expansion of the long-distance 
1—1 (LD) contributions to the form factors, but also provides important absolute normalization of the LO form factors. 
On the other hand, in the case of strong pionic decay the heavy quark remains spectator and the dynamics of the 
process is determined by the light degrees of freedom. As a consequence HQ symmetry does not provide any absolute 
normalization, thus making the pionic decays of heavy mesons a more complicated problem to analyse compared with 
, semileptonic decays. 

Q\ • The amplitudes of the strong pionic vector-to- vector and pseudoscalar-to- vector transitions have the structure 
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(V(p 2 )n(q)\P(pi)) = qvS„(p 2 )gvP7r 

(V(p 2 )n(q)\V(p 1 )) = £ v (pi)e^p 2 ) Pla p 20 e v ^g VVw , (1) 

where e denotes the polarization vector of the vector meson. The coupling constants can be expanded in a 1/mg 
series as follows 
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gvp-n = gpv-K = g ■ ' 



- 1—1 



-z-gw-K = gwix = g H V ■ . ■ , (2) 

2 THQ 

and due to the HQ symmetry in the leading 1 /toq order both constants are governed by the same quantity <?, whereas 
the higher order terms are different. 

The pionic coupling constants for heavy mesons gDD*it and gBB'ir have been analysed within various versions of 
the constituent quark model M, sum rules || and using the lattice simulations 0. 

The results of the quark model strongly depend on the particular version of the QM used: the values of g range 
from 1 (nonrelativistic quark model) to 1/3 (the Salpeter equation with massless quarks) Q. LCSR Q obtained a 
small value 0.32 ± 0.02. Recent lattice simulations reported 0.42 with however large statistical errors ±0.04 ± 0.08. 
Thus at the moment no reliable predictions are available even for the coupling constants of the ground state mesons. 
An analysis of the pionic decays of radially excited states is almost lacking. This calls for further investigation and a 
better understanding of the problem. 
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Although the present QM results are strongly scattered, one feature of this approach still seems to be very attractive: 
namely, various processes can be connected to each other via the wave functions of the mesons involved, once a proper 
formulation of the QM is used. We believe that such a proper formulation should (i) be based on a relativistic 
consideration and (ii) reproduce rigorous QCD results in the known limits, e.g. in the limit when the meson decay is 
induced by the HQ transition. The dispersion formulation of ref Q was shown to satisfy these requirements. Recently, 
we have applied this dispersion approach to analyse the B — > 7r transition form factors and determined the relevant 
quark model parameters and the wave function of the B meson Once the whole framework and the numerical 
parameters are fixed we expect to perform a reliable analysis of the pionic coupling constants of the heavy mesons. 
An additional attractive feature of this approach is that it is straight forward to incorporate the excited states, in 
particular radial excitations. 

The dispersion approach is based on the spectral representations of the form factors through the wave functions of 
the initial and final mesons. The double spectral densities are obtained from the relevant Feynman graphs, whereas 
the subtraction terms remain apriori ambiguous and should be determined from some other arguments. In ref the 
subtraction terms have been determined such that the correct structure of the 1 /toq expansion in the leading and 
subleading orders is reproduced. Such a procedure of determining the subtraction terms provides a reliable description 
of the form factors if the decaying quark is heavy for both heavy and light quarks produced in the weak decay. However 
such a strategy might be not efficient for the case when the interacting quark is light but the spectator is heavy, that 
is exactly the case of the pionic decay of a heavy meson, and requires proper modification. 

In this letter we analyse the pionic transitions of the ground state and radially excited heavy mesons, studying 
in parallel the V — » Vtt and V — > Pit cases. We make use the PCAC which allows one to reduce the calculation 
of the complicated amplitude involving three hadrons to a considerably more simple amplitudes of the P — > V and 
V — ► V transitions induced by the light-quark axial- vector current. The corresponding amplitudes have the following 
structure 

(V(p 2 )\A fl \V(pi)) = ie a ( Pl )e (p 2 )e^ a pP l/ h(q 2 )/2 + ..., 

(V(p2)|A M |P(pi)) = i£u(p2)\g„uf(q 2 ) + a+{q 2 )P ll q u + a-{q 2 )q ti q v \, (3) 

where q = pi — p 2 , P = pi + p 2 , and the dots denote terms transverse with respect to the particular form of 
which is not important for us. The form factors a_ and h contain a pole at q 2 = m 2 due to the contribution of the 
intermediate pion state in the g 2 -channel, and the residues in these poles are expressed through the pionic coupling 
constants gypu an d gvVw as follows 

h \1 ) = — T + % ) 



m„ — q* 

Vtt 

T 2 



I 2\ fir9VVTT . - / 2\ / a\ 

a -0?) = — — 2+ a -{q), (4) 



where h and a_ are regular at q 2 = to 2 . 

In the V — > V case, we find that the function $y(q 2 ) = h(q 2 )(m 2 — q 2 ) is smooth between the point q 2 = and 
q 2 = m 2 and due to the smallness of the pion mass we assume that $y(0) ~ &v(m 2 )- This yields the relation 

gvv-R = h(0)/U. (5) 
For the P — > V case an additional step is necessary: taking first the divergence of the axial-vector current we find 

$ P (q 2 ) 

< V(jp2) <fcA* \ P (Pl) >= <7i^(??2)— J 2 

m£ — q A 

with$ P (g 2 )= [f{q 2 ) + {p\ - pl)a + {q 2 ) + q 2 a^(q 2 )] (to 2 -q 2 ) + q 2 f^gvp-K- A standard PCAC assumption <E> P (to 2 ) ~ 
<!>p(0) yields 

9PV. = -^[f(0) + (pl-p 2 2)a+(0)]. (6) 

In 

Notice that the pion should not be necessarily soft and the relations (||) and @ can be readily applied to the case 
when the masses Mi and M 2 are substantially different as e.g. for the transition between the ground-state and 
radially-excited heavy mesons. 

The relations (|^) and (^) represent the coupling constants of interest through the transition form factors, for which 
we apply the dispersion approach. 



2 



I. SPECTRAL REPRESENTATIONS OF THE COUPLING CONSTANTS 



We start with considering the spectral representation of the form factors along the lines of ref [J7| . We apply this 
method to the case of the heavy spectator quark and propose proper modifications to calculate the pionic coupling 
constants. 

The dispersion approach gives the transition form factors of the meson Mi to the meson M 2 as double relativistic 
spectral representations in terms of the soft wave functions ?/>i(si) = Gi(si)/(si — M 2 ) and ^2(32) = G2(s2)/(s2— Mf ), 
of the initial and final mesons, respectively. To be specific, in the case of the weak decay Pq — > V q > induced by the 
weak quark transition Q — > q' where q' might be both heavy and light the form factor can be represented as follows 



f(l 2 ) 



dsi Gi(si) ds 2 G 2 (s 2 ) 



si 



M 2 



S2 



Ml 



f D (si,s 2 , q 2 ) + ((si - Ml) + (s 2 - Ml)) g D (si, s 2 ,q 2 ) 



+ (si - M 2 )Zi(s 1: s 2 ) + {S2 - Aff )&(*!, s 2 )] , 



(7) 



where fu 



disc Sl disc S2 /n(si, 



s 2 ,q 2 ) and gn — disc Sl disc S2 gi)(si, s 2 , q 2 ) are directly calculated from the triangle 



Feynman graph with pointlike vertices as double spectral densities of the relevant form factors. The functions £1 
and £ 2 are not known precisely but are known to behave as 1/toq in the limit toq — > 00, where tuq is the mass of 
the initial heavy quark. It is important to point out at the hierarchy of different terms in eq (H) in case the initial 



active quark is heavy: namely, the term proportional to disc Sl disc S2 /u contributes in the LO (and all other orders 
as well), the subtraction term proportional to disc Sl disc S2 3u contributes starting from the subleading order, and the 
subtraction terms proportional to £ contribute only in the higher orders. So this expansion works well if the initial 
active quark is heavy but the spectator quark is light, both for the cases of the light and heavy active quark in the 
final state. 

However in case the spectator is heavy and the active quark is light, the subtraction term proportional to §d as 
well as the functions £ contain powers of the spectator mass, and all terms have the same order of magnitude and thus 
should be considered on an equal footing. Hence, the expansion (Q) becomes ineffective in practice and the procedure 
requires a modification. 

To find such proper modification it is convenient to start with considering the V± — > V 2 transition induced by 
the axial- vector current. In the quark model the process is represented by the triangle graph with relevant vertices 
describing the quark structure of the vector meson states. The procedure described in detail in Q yields for the 
9ViV 2 tv = h(q 2 = 0)/ fir the following spectral representation 



9v 1 v 2 -!T 



mi 



OO 

- J ds^ Vl (s)^v 2 (s) 



(mi+m 3 ) 2 



A 1 / 2 1 
mi log(r) + (to 3 — mi) 1 F 

S y's + TOl + TO3 



L A V2_^i log(r) 



(8) 



Am 2 m 3 , and mi and TO3 the masses of 



with r = (s + mf - ml + A 1/2 )/(s + to? - ml - A 1/2 ), A = (s - m\ - to 2 ) 
the active and the spectator constituent quarks, respectively (we follow the notations of 0) 

This spectral representation can be rewritten in a more conventional form as an integral over the light-cone variables 
as follows 



gVxVi-n 



mi 



2tt 2 UJ x(1 



dxdk\ 



■*vi(a)*v 2 (s) 



mix + 1713(1 — x) + 



2k\ 



m 3 



(9) 



Here s = 



I "»3 



and x and 1 — x are the fraction of the light-cone momentum carried by the spectator 



X—x ' x 1 x(l- 

and the active quark, respectively. 

The expression (^) is obtained as a double dispersion representation in the invariant masses of the initial and final 
qq pairs, the spectral density of which is calculated from the triangle Feynman graph. At q 2 — it is reduced to 
a simple form (@). In principle subtraction terms can be added to this double spectral representation. However, in 
the case of the V — > V transitions there are no reasons dictating a necessity of subtractions, and we assume that the 
subtraction terms are absent. 

The representation (^) describes the ttq — > 27 decay if we set mi = m 3 and take the quark-photon vertex in the 
form q^f^q. In this case only the first term in (||) is present, and the photon wave function corresponding to the 
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pointlike interaction reads ^(s) = 1/s. So, g 117T becomes independent of the quark mass and we simply reproduce 
the value of the axial anomaly M] from the imaginary part of the triangle graph (cf |llf| and refs therein) . 

The axial- vector current satisfies the equation of motion d^A^ — 2m j'5, so the coupling constant gvti^n can be 
equivalently calculated directly from the amplitude (V2(p2)|2mj 5 (0)|Vi(pi)) = E a {p{)Ei3(p2)e f ivap<lvPvh(q 2 ) /2. 

We proceed the same way to analyse the P — * V transition: instead of calculating the full representations for / 
and a+ and then taking their linear combination @, we focus on the amplitude (V(p2)\2mj5\P(pi)). Similar to the 
V — > V case, this amplitude provides the correct double spectral density of the spectral representation for gvp-n which 
takes the form 



gvp-ir 



mi 



ds *p(s)*v(s) 



(mi+m 3 ) 2 



(s - (mi - m 3 ) 2 ) log(r) - ^1 + 



2mi 



rs + mi + m3 



((s + m 2 -m 2 )log(r)-2A 1/; 



mi 



A^UJ x{l-xf 



dXdkl *p(*)*v(») 



s - (mi - m 3 ) 2 - 



l-x 



1 



2m 1 



's + mi + m3 



(10) 



The form of the relations (^|) and (Q) however prompts that in distinction to the V — > V transition, in the P — > V 
transition the subtraction term is nonzero. The latter cannot be determined uniquely within the dispersion approach. 
It is very important however that the HQ symmetry ensures the subtraction term to contribute only in the subleading 
1/mQ order. Namely, the HQ symmetry predicts the LO relation between the coupling constants gvp-ir = n^Qgw-K- 
The double spectral densities of the representations for gvv-?r ® and gvPn (0) satisfy this relation. Hence the 
subtraction terms in gvv-n and gyp^ should also have the same LO 1/mQ behavior. Since the subtraction term in 
gvv-K is absent, the subtraction term in gvp-n does not contribute in the LO. Although we cannot determine the 
subtraction term uniquely, several reasonable ways of fixing this subtraction term yield a numerical uncertainty in 
g PVrr to be not more than 10%. 

The normalization condition for the radial wave functions ^ of the ground state and the radial excitation of the 
vector and the pseudoscalar mesons has the form 



1 f A 1 / 2 

ds Ms) — - 

57T Z J S 



(mi - m 3 y 



1 

8^ 



dxdk 2 



x(l 



(mi - m 3 ) 2 ] = S k 



(11) 



It is easy to see that in the nonrelativistic (NR) limit \k\ = A 1 / 2 /2y / s <C m\,m§ the coupling constants take the 
values gpvK = 2A-/// ff and gyy^ — 2/ f n . Moreover, in the NR limit the coupling constants of the transition between 
the ground state and the radial excitation vanish, i.e. gpyi^ — gyv'-n = 0- 

For the analysis of the HQ expansion it is convenient to introduce a new variable z such that s = (z + mi + m 3 ) 2 
and use the fact that the soft wave functions </>o are localized in the region z ~ Aqc d ■ Performing the HQ expansion 
of all quantities including 4" in the inverse powers of m 3 = ttiq — ► 00 and keeping mi = O(Aqcd) (see 0) we find 
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9ij 



dz(f)l(z)4f' Q (z)mi 



, z + mi + \/ z(z + 2mi ) \ , — ; 

mi log I -^-^ V ) = U ) + 2^z(z + 2 mi ) 



z + mi — y z(z + 2m\) 



(12) 



In this expression Vo 
the normalization condition 



J are the LO soft wave function of the ground state and the radially excited mesons which satisfy 



dzcf>i(z)4i(z)y/z'(z + 2mi) 3/2 = S i: 



(13) 



Notice that the Isgur-Wise function is directly expressed through <j> (see 0]). These formulas provide a possibility to 
evaluate the coupling constants of interest if the numerical parameters of the model are known. 



II. NUMERICAL ANALYSIS 



We now provide the numerical estimates of the LO quantity g for the ground-state and radially excited mesons. The 
spectral representation for this quantity is completely fixed and we need to choose the proper numerical parameters 
of the model. As found in many applications of the dispersion approach (see and refs therein), an approximation 
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of the soft wave function by an exponent provides reasonable estimates. Assuming the wave function to have the form 
^(s) ~ exp(— fc 2 /2/3 2 ), the LO radial wave function of the ground state reads Q 



(z) "^7T^ exp l 2^—'- (14) 



Similarly, we assume for the wave function of the radial excitation the form 

Mz) * I 1 " Cr ^-) exp I — wT) ' (15) 

where the normalization factor and the coefficient C r are fixed by the normalization conditions and the orthogonality 
between the wave functions of the ground state and the radial excitation. 

The light-quark mass was determined from the numerical matching of the transition form factors in the B — > tt. p 
decay to the available lattice data (|] and was found to be tightly restricted to the range m\ — 0.23 ± 0.01 GeV '. 
Therefore the only unknown parameter is f3 oo- The analysis of the leptonic decay constants of heavy mesons shows 
that simple approximate relations 

0p{mo) = /3oo(l - Cp/mo), /V( m Q) = _ CV/ m Q) ( 16 ) 

with Cp ~ 0.1, Cv — 0.2 and /?oo = 0.5 yield the values of (3 which provide reasonable leptonic constants of the 
charm and beauty mesons (see Table 1) calculated within the quark model through formulas given in 0,|l2|. The 
slope parameter of the IW function for such (3^ is p 2 = 1.2 ± 0.03 0. The eq ( [l2| ) then yields 

g = 0.5, g r = 0.11. (17) 

Thus for the LO quantity the orthogonality of the radial wave functions provides a strong suppression of g r ^. 

The results of calculating the pionic coupling constants of the B and D mesons are given in Table 2. Whereas the 
gvvir is determined quite reliably within the dispersion apprach, the gpvn suffers from the intrinsic uncertainty of the 
approach based on the impossibility to fix the subtraction terms. Assuming several reasonable choices of subtraction 
terms in the form factors /, a+ and g which provide the correct behaviour in the limit of the heavy active quark, 
but differ for the case of the heavy spectator we have found that in all cases, the difference in the coupling constants 
of the ground state is not more than 10%. This minor difference is easily understood taking into account that the 
subtraction terms contribute only in the subleading order. 

Notice that there is a possibility of an alternative determination of the gp* p-n from the form factors of the semilep- 
tonic B — > 7r transition near zero recoil. Namely, 

f (M 2 B ) ~ / (A/* 2 ) = ggg^ (1 + r'(l)), (18) 



where r(q 2 = q 2 /M^,) is defined as follows 



f-(q 2 ) m 



2 



The function r(q 2 ) satisfies the condition r(l) = 1 and is known numerically from the results of the lattice simulations 
in the region q 2 < 0.7 0, in particular r(0.7) ~ 0.9. The function r{q 2 ) is regular near q 2 = 1 and thus there are no 
reasons for fast variations of this function near q 2 = 1 . Then, taking into account the lattice data we can estimate 
r'(l) < 0.8. Using the current algebra relation /o(M^) = fp/fn, neglecting the 1/mt, corrections which are small, 
and taking the ratio of the leptonic constants fp* / fp ~ 1.2 wc find 

1 > 0.45 (20) 



/W/s(l + r'(l)) 



1 Notice that in contrast to the PCAC based approaches, the NR quark- model calculations based on taking into account the 
quark structure of the emitted pion through the pion wave function yield bigger values of g and do not yield a suppression of g r 
since the product of the orthogonal wave function is smeared by the integration with the pion wave function and as the result 
the orthogonality is not efficient p3|. 
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in agreement with ( J17| ) and with the recent result of the lattice simulation Jd| . We would like to point out that the 
value g ~ 0.27 proposed in jlj] as the preferrable solution found from the description of the D* decays falls far out 
of this estimate. 

The calculation of the coupling constants of the pionic transition between the ground state and the radial excitation 
is more involved. In this case the LO term is itself strongly suppressed because of the orthogonality of the wave 
functions. So the dependence on the parameters of the wave function of the radial excitation in the case of the 
gvv'TT and, in addition to this, also the dependence on the particular choice of the subtraction procedure in gpv'ir is 
more sizeable. Table 2 presents numerical results. The slope parameter of the Gaussian wave function of the radial 
excitation at finite masses is not known, so for estimating gvv-n we assume that this parameter lies within a 10% 
interval around the corresponding ground-state slope. In the case of gpv'-n the relation @ prompts that the 1/toq 
corrections numerically might be sizeable if the M v , — Mp is big as it is for the radial excitations of heavy mesons. 
In this case we cannot obtain a reliable estimate and Table 2 provides the lower bounds found from the spectral 
representations without subtractions. 



III. CONCLUSION 



We have analysed the pionic coupling constants of heavy mesons within the framework based on the combination 
of PCAC with the dispersion approach and obtained the following results: 

1. Spectral representations of the coupling constants in terms of the wave functions of the initial and final heavy 
mesons have been obtained. The double spectral densities of these spectral representations for gpv-n and gvv-n are 
equal in the HQ limit in agreement with the HQ symmetry. 

In the case of the V — ► V transition the calculation of the gvv-ir is equivalent to the calculation of the amplitudes 
(V|2mj5|V). There is no reason dictating a necessity of any subtraction term in the spectral representation for gvvn 
and thus the latter is determined unambiguously within our approach. 

In the P — > V case the double spectral density of gvPn can be also found from the amplitude (V\2mj5\P). The 
spectral representation for gvp-n, however, contains a subtraction term which cannot be fixed within our approach. 
Important is that due to the HQ symmetry this subtraction term does not contribute in the leading 1 /itiq orderQ and 
thus numerically is not essential at least for the gs* Bn- 

2. A spectral representation for the LO 1/toq coupling constant g which describes the pionic transition in the HQ 
limit has been obtained. The g is represented through the LO wave functions of the heavy mesons which determine 
also the Isgur-Wise function. The g is determined quite reliably within the dispersion approach since it is not affected 
by the uncertainties in the subtraction procedure. 

For the transition between the ground-state heavy mesons we have found g = 0.5±0.02. We argue that the behavior 
of the B — > 7r transition form factors near zero recoil provide independent arguments in favour of such considerably 
large value of g. 

For the pionic transition between the ground state and the radial excitation in the HQ limit a strong suppression 
g r = 0.11 ± 0.02 is observed as a consequence of the orthogonality of the corresponding radial wave functions. 

3. Using the QM parameters obtained from the analysis of the B — > it decay we have calculated the pionic coupling 
constants of the ground state charm and bottom mesons. Our final numerical estimates are listed in Table 2. 

For the V — > V transition higher-order I/tuq effects in gvv-n are found to be small and not to exceed 5-6%. 

In the P — > V case, the coupling constant gpvn depends on the choice of the subtraction procedure which cannot 
be fixed unambiguously. However, numerically the uncertainty estimated by using several reasonable subtraction 
procedures is not more than 10%. 

4. We estimated the pionic couplings of the radially excited states B* and D* . The errors turn out to be much 
bigger compared with the corresponding ground states: the LO contribution is strongly suppressed and thus the 1/toq 
effects are found to be numerically more important. As a result the coupling constants are sensitive to the specific 
form of the radial-excitation wave function and gpv'-n, in addition to this, strongly depends on the the details of the 
subtraction procedure. Therefore only lower bounds on gpv'-n are given. 



2 Strictly speaking, we cannot completely exclude a possibility of subtraction terms both in jyv, and gpvrr which give the 
same contributions to both quantities in the LO. However, unless no reasons for such subtraction terms are found we do not 
include them into consideration. In a recent analysis of the B — > ty form factors [^| we have found that the behavior of fo at 
zero recoil is compatible with gs* b» — 0.5 — 0.7. The values of gs* Btt as big as 0.6 — 0.7 can be obtained only by assuming the 
subtraction terms both in gB*B-n and gs*B*-n which give a nonvanishing contribution already in the leading l/mq order. 
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TABLE I. Quark masses and the slope parameters of the soft meson wave functions and the calculated leptonic decay 
constants (in GeV) 





m c 


m u 


/3oo 


Pb 


Pb* 


(3d 


Pd* 


fB 




U 




4.85 


1.4 


0.23 


0.5 


0.49 


0.49 


0.46 


0.43 


0.16 


0.175 


0.2 


0.22 



TABLE II. Pionic coupling constants of heavy mesons. The error bars correspond to the variations of the slope parameters 
f3 around the average values of Table 1 yielding a 10% variations of the leptonic decay constants. 















gvvir 


9vptv 


gvviT 


gvpjr 


D,D* 


0.53L0.03 


0.53L0.05 


0.15L0.03 


> 0.14 


B, B* 


0.5L0.02 


0.5L0.04 


0.12±0.03 


> 0.12 


HQ-limit 


0.5±0.02 




0.11L0.02 
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